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Abstract -Satellite remote sensing estimates of Gross Primary Production (GPP) have routinely 
been made using spectral Vegetation Indices (Vis) over the past two decades. The Normalized 
Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), the green band 
Wide Dynamic Range Vegetation Index (WDRVI gre en), and the green band Chlorophyll Index 
(CIgreen) have been employed to estimate GPP under the assumption that GPP is proportional to 
the product of VI and photosynthetically active radiation (PAR) (where VI is one of four Vis: 
NDVI, EVI, WDRVIgreen, or CIgreen)- However, the empirical regressions between VI*PAR and 
GPP measured locally at flux towers do not pass through the origin (i.e., the zero X-Y value for 
regressions). Therefore they are somewhat difficult to interpret and apply. This study 
investigates (1) what are the scaling factors and offsets (i.e., regression slopes and intercepts) 
between the fraction of PAR absorbed by chlorophyll of a canopy (fAPAR c hi) and the Vis, and (2) 
whether the scaled Vis developed in (1) can eliminate the deficiency and improve the accuracy 
of GPP estimates. Three AmeriFlux maize and soybean fields were selected for this study, two of 
which are irrigated and one is rainfed. The four Vis and fAPAR c hi of the fields were computed 
with the MODerate resolution Imaging Spectroradiometer (MODIS) satellite images. The GPP 
estimation performance for the scaled Vis was compared to results obtained with the original Vis 
and evaluated with standard statistics: the coefficient of determination (R ), the root mean 
square error (RMSE), and the coefficient of variation (CV). Overall, the scaled EVI obtained the 
best performance. The perfonnance of the scaled NDVI, EVI and WDRVIgreen was improved 
across sites, crop types and soil/background wetness conditions. The scaled CI gre en did not 
improve results, compared to the original CIgreen- The scaled green band indices (WDRVIgreen, 
CIgreen) did not exhibit superior performance to either the scaled EVI or NDVI in estimating crop 
daily GPP at these agricultural fields. The scaled Vis are more physiologically meaningful than 
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original un-scaled Vis, but scaling factors and offsets may vary across crop types and surface 
conditions. 

Key Words - Daily GPP, MODIS, Vegetation Index, fAPAR C hi 


I. INTRODUCTION 

Atmospheric general circulation models require quantification of land-atmosphere 
exchanges of energy, water and momentum, including CO 2 fluxes which can be provided by land 
surface process models (Bonan et ah, 2011; Dickinson et ah, 1993; Sellers et ah, 1986). Satellite 
remote sensing offers inputs such as land cover types and the Nonnalized Difference Vegetation 
Index (NDVI) (Deering, 1978; Tucker, 1979) for use in the land surface modeling (Dickinson et 
al., 1990; Sellers et al., 1994). Pioneering work (Asrar et al., 1992; Myneni et al., 1997; Running 
et al., 2000; Sellers, 1987) has shown the fraction of photosynthetically active radiation (PAR) 
absorbed by a canopy/vegetation (FPAR, i.e., fAPAR canop y) can be approximated with NDVI 
(Running et al., 2000). Therefore, NDVI has been employed to estimate vegetation Gross 
Primary Productivity (GPP) in a variation (as G P P=e* NDVI* PAR, Running et al., 2000), 
inspired by the logic from the Light Use Efficiency (LUE) model (Monteith, 1972, 1977): 

GPP = s * fAPARpsN * PAR = s *APARp S N, (1) 

where e is LUE for vegetation photosynthesis (PSN) (Running et al., 2000) and fAPARpsN is the 
fraction of PAR absorbed for PSN (APARpsn)- Monitoring changes in crop GPP with satellite 
remote sensing data advances the capability to understand and manage global food security, 
sustainability practices, and environmental impacts, and to study global carbon cycle and global 
water cycle. 
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The three-band Enhanced Vegetation Index (EVI) (Huete et ah, 1997) and the two-band 
EVI (called EVI2, Jiang et ah, 2008) have also been utilized to predict terrestrial GPP in a 
similar way as GPP=£*EVI*PAR (Jin et ah, 2013; Kalfas et ah, 2011; King et ah, 201 1; Li et ah, 
2007; Mahadevan et ah, 2008; Schubert et ah, 2012; Sjostrom et ah, 2011; Wu et ah, 2008, 2010, 
2011, 2012; Xiao et ah, 2004; Yan et ah, 2009). In addition, Gitelson and colleagues also 
explored the application of the green band Wide Dynamic Range Vegetation Index (WDRVIgreen) 
and the green band Chlorophyll Index (CIgreen) for crop GPP estimation, in addition to the NDVI 
and EVI (Gitelson et ah, 2008, 2012; Peng and Gitelson, 2011, 2012; Peng et ah, 2011). 

However, since the empirical regressions between the VI*PAR products and GPP 
measured locally at flux towers do not pass through the origin (i.e., the zero X-Y value for 
regressions) and produce offsets, they are somewhat difficult to interpret and apply (Gitelson et 
ah, 2012; Sims et ah, 2006; Zhang et ah, 2014b). This is considered to be a source of error 
affecting the accuracy and reliability of remote sensing GPP estimates based on Vis. In the 
literature, there is no paper that presents how to scale the Vis in space and time to solve the 
problem. 

The standard MODerate resolution Imaging Spectroradiometer (MODIS) 8-day GPP 
product (MOD17A2 GPP) uses the MOD15A2 FPAR (a fAPAR canopy ) product as a model input 
(Running et ah, 2004; Zhao and Running, 2008). Investigations to find the scaling factor and 
offset of NDVI through fAPAR canopy - NDVI functions have been conducted, where fAPAR canopy 
=a 0 *NDVI+b 0 (a 0 is the scaling factor or slope, and bo is y intercept or offset) (Fensholt et ah, 
2004; Goward and Huemmrich, 1992; Knyazikhin et ah, 1998, 2002; Potter et ah, 1993; Prince 
and Goward, 1995; Randerson et ah, 1996; Sellers et ah, 1996; Sims et ah, 2005). However, the 
MOD15A2 FPAR product overestimates in-situ fAPAR canopy during spring greenup and fall 
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senescent periods, and underestimates in-situ fAPAR can opy in mid-summer during peak GPP 
activity at the agricultural fields we selected [see (Zhang et ah, 2014a) for details]. 

We developed an algorithm to retrieve the fraction of PAR absorbed by chlorophyll 
throughout the canopy (fAPAR c hi) from actual MODIS observations or from synthesized 30 m 
MODIS-spectral-like observations simulated with EO-1 Hyperion images (Zhang, 2003; Zhang 
et ah, 2005, 2009, 2012, 2013,2014c). We found that fAPAR cM ^ fAPAR canopy , and that the 
fraction of PAR absorbed by foliage non-chlorophyll components (fAPAR non _ C hi) varies with 
types and seasonally (Zhang et ah, 2013). Zhang et al. (2014a) presented the performance of 
£APAR c hi and MOD15A2 FPAR in crop GPP estimation, and concluded that fAPARchi is 
superior to MOD15A2 FPAR. Zhang et al. (2014b) investigated the performance of original un- 
sealed Vis in GPP estimation, and suggested that further investigation on the performance of 
scaled Vis should be carried out. 

The objectives of this paper are straightforward: 1] to explore how surface conditions 
affect the scaling factors (“a”) and offsets (“b”) derived through regression analysis of fAPAR c hi 
vs. the four Vis: fAPAR c hi =a*VI+b for each crop type (com, soybean) per field; 2] to investigate 
how much the scaled Vis can improve the prediction accuracy of GPP estimates compared to the 
prediction of original un-scaled Vis. 


II. METHODS 

II. 1 Study sites and tower data 

The three AmeriFlux crop sites for com, or maize (Zea mays L.) and soybean ( Glycine 
max [L.] Merr.) used in this study are located at the University of Nebraska-Lincoln (UNL) 
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Agricultural Research and Development Center near Mead, Nebraska (US-NE1, US-NE2 and 
US-NE3). The first two fields are circular (radius ~ 390 m) and equipped with center-pivot 

irrigation systems (US-NE1, 41°09’54.2”N, 96°28’35.9”W; US-NE2, 41°09’53.6”N, 
96°28’07.5”W). The third is a 790 m long square field (US-NE3, 41°10’46.7”N, 

96°26’22.4”W) that relies entirely on rainfall. Each field is equipped with an eddy covariance 
flux tower (Gitelson et ah, 2012; Gitelson et ah, 2006; Peng et ah, 2013). The first field (US-NE1) 
is a continuous maize field while the other two fields are maize-soybean rotation fields (soybean 
is planted in even years). 

Tower eddy-covariance carbon exchange, PAR, and GPP measurements in growing 
season from 2001- 2006 are publically available and can be downloaded from 
ftp://cdiac.oml.gov/pub/ameriflux/data. The nighttime ecosystem respiration/temperature Qio 
relationship was used to estimate the daytime ecosystem respiration (Baldocchi, 2003). Daily 
GPP was computed by subtracting respiration (R) from net ecosystem exchange (NEE), i.e., 
GPP=NEE-R (Suyker et al., 2005). These sites provided the opportunity to examine the semi- 
empirical relationships between fAPARchi versus Vis for both C4 (maize) and C3 (soybean) 
crops in both irrigated and non-irrigated ecosystems, and to investigate the benefits of employing 
the scaled relationships to estimate GPP. 

II. 2 Remote sensing data processing and GPP estimation 

Six years (2001-2006) of MODIS LIB calibrated radiance data (MOD021KM and 
MOD02HKM) and geolocation data (MOD03) covering the three study sites were downloaded 
from https://ladsweb.nascom.nasa. gov:9400/data/ . Two of the MODIS bands have a nadir spatial 
resolution of 250 m: B1 (red, 620 - 670 nm) and B2 (near infrared, NIRi, 841 - 876 mn). The 
MODIS land bands 3-7 have a nadir spatial resolution of 500 m: B3 (blue, 459 - 479 nm), B4 
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(green, 545 - 565 nm), B5 (NIR 2 , 1230 - 1250 nm), B6 (shortwave infrared, SWIRj, 1628 - 
1652 nm) and B7 (SWIR?, 2105 - 2155 nm). The centers of the original 500 m grids defined in 
the standard surface reflectance product (MOD09) that encompass the three tower sites are not 
the centers of the three fields and vegetation in each of the original 500 111 grids is not 
homogeneous [see Figure 2 of (Guindin-Garcia et ah, 2012)]. The MODIS gridding procedure 
for the standard MOD09 product does not ensure the gridded surface reflectance covers the 
entire grid (Wolfe et ah, 1998). A modified gridding procedure was used for this study (Zhang et 
al., 2014b), whereby the centers of the three 500 m grids were matched to the centers of the three 
fields, respectively. The LIB radiance data from each swath were gridded at 500 m resolution for 
MODIS bands 1-7 with area weight of each MODIS observation. This modified gridding 
processing was incorporated into the Multi- Angle Implementation of Atmospheric Correction 
(MAIAC) algorithm (Lyapustin et al., 2008, 201 la, b, 2012). MAIAC is an advanced algorithm 
which uses time series analysis and a combination of pixel-based and image -based processing to 
improve cloud/snow detection, and to achieve more accurate aerosol retrievals and atmospheric 
correction, based on the bidirectional reflectance distribution function (BRDF) model of the 
surface. 

Derived bidirectional reflectance factors (BRF, also called directional surface reflectance) 
in MODIS bands 1-7 were used for this study. The impact of MODIS observation footprint size 
resulting from variable view zenith angle (VZA) on crop daily GPP estimation for these sites 
was recently reported elsewhere (Zhang et al., 2014b). In order to eliminate the potential bias due 
to large VZAs, only observations with VZA < 35° were included in this study. The surface 
reflectance data (p) were used to calculate the following indices (Deering, 1978; Gitelson, 2004; 
Gitelson et al., 2007, 2012; Huete et al., 1997, 2002; Tucker, 1979): 
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We used the PROSAIL2 model (Jacquemoud and Baret, 1990; Baret and Fourty, 1997; 
Braswell et al., 1996; Verhoef, 1984, 1985; Zhang et ah, 2005, 2009, 2012, 2013), a coupled 
soil-canopy-leaf radiative transfer model, to retrieve fAPARchi, the fraction of PAR absorbed by 
the foliage of the canopy (fAPARf 0 ij age ), and the fraction of PAR absorbed by the non- 
photosynthetic foliage components (fAPAR non _ cM )(Zhang et ah, 2014a). A pixel is composed of 
canopy and soil (Zhang et al, 2009, 2012, 2013). The canopy is partitioned into foliage and stem 
(including branch), and the foliage component is further partitioned into chlorophyll (chi) and 
non-chlorophyll (non-chl) components, where non-chl is composed of non-photosynthetic 
pigments (referred to as brown pigment) and dry matter (Baret and Fourty, 1997). The surface 
reflectances of MODIS bands 1 - 7 are used for retrieval of 1APAR variables (Zhang et al., 2009, 
2012, 2013,2014c): 

fAPAR non-chl = fAPAR brownpigment + fAPAR dry _ matter (6) 

fAPAR foliage = fAPAR chl + fAPAR non _ chl (7) 

fAPAR canopy = fAPAR foliage + fAPAR stem (8) 
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The scaling factors (“a”) and offsets (“b”) of Vis were derived from linear regression through 
fAPAR c hi - VI functions for each crop type per field, where fAPAR c hi =a*VI+b (VIs=NDVI, 
EVI, WDRVIgreen, and CIgreen). 

The product of Vis and tower daily PAR (VI*PAR) and the product of scaled Vis and 
daily PAR (scaled VI* PAR) were compared against the tower daily GPP for each crop type per 
field (GPP=^*VI*PAR or GPP=f ^scaled VI*PAR). The coefficients and “e ” were 
computed with a least squares best fit algorithm. The computed values for ^ and £ were then 
used to predict GPP, and coefficient of determination (R~), the root mean square error (RMSE, g 
C m'“ d' ) and coefficient of variation (CV, %) was calculated. The average light use efficiency 
at chlorophyll level (LUE C hi, i.e., £ cht ) was computed using GPP=LUE c hi*fAPAR c hi*PAR with a 
least squares best fit algorithm. Improvements of crop daily GPP estimation using scaled Vis 
were assessed. 


III. RESULTS 

The scaling factor (“a”, also called slope) and offset (“b”, also called y-intercept) 
obtained through the regression functions fAPAR c hi =a*VI+b for each crop per field are listed in 
Table 1, where the statistics for the R , RMSE and x-intercept are also summarized. The x- 
intercepts of fAPAR c hi =a*VI+b give minimum VI values at zero fAPARchi. The 95% confidence 
intervals of slope, y-intercept and x-intercept for each crop per field are reported, too. The CIgreen 
is a simple ratio index while the other three Vis include consideration of normalization. The 
confidence intervals for CI gre en are different from those for other three Vis for each type per field. 
For each crop type in irrigated fields USNE1 and USNE2, the confidence intervals of y- 
intercepts and x-intercepts for NDVI, EVI and CIgreen are different from each other. For each 
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crop type in rainfed field USNE3, the confidence intervals of y-intercepts and x-intercepts for 
NDVI and CI gre en overlap each other, but are different from those for EVI. Mean values of the 
confidence intervals of the slopes, y-intercepts and x-intercepts vary with Vis, sites, crop types 
and irrigation options. None of the y-intercepts or x-intercepts for NDVI, EVI or WDRVIgreen is 
close to the origin (i.e., zero X-Y point). 

The functions in Tab. 1 were used to compute the scaled values of NDVI, EVI, 
WDRVIgreen and CIgreen for each crop type per field. For instance, for the NDVI at US-NE1 : 
scaled NDVI = 1.1 l*NDVI-0.29. The coefficients"^ and e and LUE c hi of each crop per field are 
listed in Table 2. Corn LUE c hi is ~1.6 times of soybean LUE c hi(Tab. 2), which agrees with the 
expectation that C4 plants have higher LUE than C3 plants (e.g., Prince, 1991), and explains why 
maize displays a wider daily GPP range (~34 g C m" d" ) than soybean ( (-19 g C m" d' )(Zhang 
et ah, 2014b). The coefficients’^ and e were applied to estimate crop daily GPP. 

Figure 1 shows the estimated soybean daily GPP for the rainfed field US-NE3 using the 
four original Vis with T^ and the scaled Vis withT, compared to tower daily GPP. The scaled 
NDVI, EVI and WDRVI gree n combined with T had better GPP estimation perfonnance than the 
original counterparts, respectively, demonstrating higher R and lower RMSE. Compared to the 
original counterparts, the (scaled NDVI)*PAR, the (scaled EVI)*PAR and the (scaled 
WDRVIgreen) 515 PAR values were closer to 0 when GPP=0. The scaled CI gree n did not provide 
better GPP estimation than the original CIgreen- In order to save pages, similar figures for US-NE1, 
US-NE2 and figures for maize in US-NE3 are not presented in this paper. 

Table 3 summarized the statistics (R , RMSE and CV) for estimating crop daily GPP 
using the original Vis with ~e q and the scaled Vis with £, respectively. These statistics show that 
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the best performance was obtained with the scaled EVI while the least successful performance 
among the four scaled Vis was obtained with CI gre en across the sites, crop types and irrigation/ 
rainfed options. For example at the US-NE1 site, scaled EVI and scaled CI gree n had contrasting 
best/worst performances in GPP estimation: R 2 : 0.88/ 0.77, RMSE: 2.92/4.05 g C m' 2 d" 1 , and 
CV: 19%/ 26% (Tab. 3). GPP estimates for corn had better perfonnance than for soybean using 
scaled NDVI and EVI for sites US-NE2 and US-NE3. Better results might be achieved for the 
sites examined in other studies (King et ah, 2011; Sjostrom et ah, 2009) if the scaled EVI 
(through coefficients obtained from the regression of fAPARchi vs. EVI) had been utilized. 

For each crop in any field, the scaled NDVI, EVI and WDRVIgreen improved the 
prediction performance of crop daily GPP while the scaled CIgreen did not, compared to the 
original un-scaled Vis. GPP improvements for the three that benefited from scaling, ranked from 
most to least were the NDVI, WDRVIgreen, EVI, for which the R 2 increased (f: 0.16, 0.13, 0.09), 
RMSE decreased (f:0.95, 0.78, 0.65 g C m' 2 d' 1 ), and the CV also decreased (|:8%, 6%, 5%). 

The improvements also varied with crop types and irrigation conditions. For example, the NDVI 
improvement for soybean (R 2 , ) 0.20; CV, | 9%) was better than for corn (R 2 , |0.13; CV,|7%), 
and the average improvement for the rainfed field (R 2 ,|0.2 1 ; RMSE, jl.10 g C m' 2 d" 1 ; and CV, 

1 10% ) was better than for the irrigation fields (R 2 , t0. 12; RMSE, 4,0.85 g C m' 2 d" 1 ; and 
CV,|6%). 


IV. DISCUSSION 

The PSROAIL2 model well distinguishes vegetation from soil and fAPAR c hi retrieved 
with the PROSAIL2 model excludes the impact of soil/background (Zhang et ah, 2012, 2013). 


11 


252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 


The fAPARf 0 ii a ge comprises chlorophyll and non-chlorophyll foliage fractions (fAPAR c hi, 
fAPAR n on-chi). Therefore, the PAR absorbed by non-photosynthetic vegetation components (NPV) 
of the canopy is excluded from APAR c hi since APAR c hi=fAPAR c hi*PAR This is the theoretical 
basis for potential improvement of GPP estimation using the scaled Vis. The x-intercept values 
of the semi-empirical linear functions of fAPAR c hi vs. VI in Table 1 have an important 
biophysical meaning: there is not any chlorophyll showing up at the pixel when its un-scaled VI 
is less than its x-intercept value. Gitelson and colleagues (Gitelson et ah, 2007) reported that, 
before green-up when green leaves do not appear, MODIS 250 m NDVI values for the fields 
could be greater than 0.2, which is close to the minimum x-intercepts of NDVI (0.23, Tab. 1) we 
found with MODIS 500 m images. In irrigated fields, the mean values of the x-intercept 
confidence intervals for EVI were about half of those for NDVI, and about 1/3 as large as those 
for WDRVIgreen (Tab. 1). In rainfed fields, the mean values of the x-intercept confidence 
intervals for EVI were about half of those for both NDVI and WDRVI gree n (Tab. 1). 
Soil/background wetness has less impact on EVI than on NDVI which is consistent with the 
original idea that inspired the development of EVI (Huete, 1988; Huete et ah, 1997). Daughtry et 
al. (2000) has expressed that Vis combined with NIR and red bands are less impacted by 
background than Vis combined with NIR and green bands. Earlier studies (Sims et ah, 2006, 
2008) have shown that GPP drops to zero at variable EVI values (i.e., x-intercept EVI values) in 
their selected flux sites, and have found the minimum x-intercept value is ~0.1. So Sims et al. 
(2008) has developed a GPP model using EVI -0.1 instead of the original EVI. The x-intercept 
confidence intervals of EVI in the three fields (US-NE1, US-NE2 and US-NE3) ranged from 
(0.12, 0.13), (0.14, 0.15) to (0.16, 0.18). Our findings are consistent with earlier empirical studies 
(Daughtry et al., 2000; Huete, 1988; Huete et al., 1997; Sims et al., 2006, 2008). Furthermore, 
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the scaled Vis with scaling factors and offsets using the semi-empirical relationships between 
fAPAR c hi vs. Vis for each crop type per field are more physiologically meaningful (Tab. 1) than 
the original un-scaled Vis. 

The ~£ estimates for all scaled Vis are close to the relevant LUE c hi values for each crop 
type per field. In contrast, the estimates associated with the original un-scaled NDVI and 
WDRVIg reen are lower than the relevant LUE c hi values. The estimates for CIg reen are much 
lower than the relevant LUE ch i values because the original un-scaled CI gre en range (~ 1 to 10) is 
much wider than the scaled CIgreen range (~0 to~l). It is worth noting that both the ~e^ and the T 
estimates for the original EVI and the scaled EVI are close to the physiologically relevant LUE c hi 
values. This partly explains the reasonableness and success of the Vegetation Photosynthesis 
Model (VPM) (Xiao et ah, 2004) which assumes GPP=e*EVI*PAR. This study suggests that the 
GPP estimation made with the VPM may be improved by replacing the original EVI with 
fAPAR c hi, or by scaling the EVI using the relationship between fAPAR c hi and EVI. 

The R between tower daily GPP and estimated GPP with scaled Vis for all cases ranges 
from 0.66 to 0.88 while the RMSE (CV) between them ranges from 4.37 to 2.1 1 g C m' 2 d' 1 
(from 31% to 17%). Although the R 2 between fAPAR ch i and scaled VI is high for all cases (0.73 
- 0.97), the RMSE between fAPAR ch i and scaled VI varies with crop type, irrigation/rainfed 
options, and VI options, which caused the variation of the performance of estimated GPP with 
scaled Vis. Among the four scaled Vis, the RMSE between fAPARchi and the scaled EVI is 
smallest and the R 2 is highest for all study sites. For US-NE2 and US-NE3, the RMSE between 
fAPARchi and scaled CIgreen is biggest and the R" is lowest. 
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V. CONCLUSION 


This study exhibited improvement in the performance of crop daily GPP estimation using 
scaled NDVI, EVI and WDRVIgreen, compared to their original un-scaled counterparts. However, 
performance improvement of crop daily GPP estimation using scaled CI gre en was not observed. 
The irrigated fields have better performance, as compared to the rainfed field. The performance 
also varied with crop types and VI options. The scaled EVI provided the best performance 
among all cases. This study does not find that the scaled WDRVIgreen or the scaled CIgreen is 
superior to the scaled NDVI or scaled EVI in predicting crop daily GPP. 

Compared to the original Vis, the scaled Vis developed with the semi-empirical 
relationships between fAPAR c hi and Vis are more physiologically meaningful. However, the 
scaling factors and offsets (and x-intercepts) vary field by field, and vary type by type. 
Investigations to explore the scaling factors and offsets of these Vis using fAPAR c hi for other 
plant functional types should be carried out in the future. We will explore how the scaling factors 
and offsets change over space and time, and vary with climate. Investigations on whether scaled 
EVI is best for all fields and all types among the four scaled Vis are also needed. We suggest an 
approach whereby MODIS-derived Vis are scaled pixel by pixel. This approach provides scaled 
Vis for use when fAPAR c hi is unavailable. We expect that future research on GPP simulation 
based on the biochemical or land surface modeling (Bounoua et al., 2000; Potter et al., 2003; 
Sellers et al., 1994, 1996) will achieve reduced uncertainty and improved accuracy when the 
scaled MODIS Vis replace the original Vis. 
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Figure captions 

Figure 1. Comparison between tower daily GPP vs. estimated daily GPP for the US-NE3 site 
(soybean): (a) NDVI; (b) EVI; (c) WDRVI green ; and (d) CI gree n- Filled circles use original un- 
sealed Vis while empty circles use scaled Vis. Only observations with VZA < 35° are included. 

Table captions 

Table 1. List of relationships between fAPAR c hi and Vis for the three crop sites (y=ax+b, 
y:fAPAR ch i, x:VI). The 95% confidence intervals of slope (“a”), y-intercept (“b”), and x- 
intercept are presented. Coefficients of determination (R ) and root mean square error (RMSE) 
are also presented. 

Table 2. List of the coefficient e~ 0 in GPP=£ji*VI*PAR, the coefficient e in GPP=£ ^scaled 
VI* PAR, and LUE ch , in GPP=LUE ch i*fAPAR ch i*PAR (unit: g C mol' 1 PPFD) 

Table 3. Coefficients of determination (R ), root mean square errors (RMSE, g C m' d' ) and 
coefficients of variation (CV) for simulated GPP with the Vis using two options: original 
unsealed Vis versus scaled Vis, compared to tower daily GPP 
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Estimated GPP (g C rrr 2 d 1 ) Estimated GPP (g C rrr 2 d 1 ) 

Figure 1. Comparison between tower daily GPP vs. estimated daily GPP for the US-NE3 site 
(soybean): (a) NDVI; (b) EVI; (c) WDRVI gree n; and (d) CI gree n. Filled circles use original un- 
sealed Vis while empty circles use scaled Vis. Only observations with VZA < 35° are included. 
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Table 1. List of relationships between fAPAR c hi and Vis for the three crop sites (y=ax+b, y:fAPAR C hi, x:VI). The 95% 
confidence intervals of slope (“a”), y-intercept (“b”), and x-intercept are presented. Coefficients of determination (R~) and 
root mean square error (RMSE) are also presented. 
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Table 2. List of the coefficient £ 0 in GPP=£ 0 *VI*PAR, the coefficient £ in GPP=£ ^scaled VI*PAR, and LUE c hi in 
GPP=LUE ch i * fAP ARchi * PAR (unit: g C mol' 1 PPFD) 
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2 2 1 

Table 3. Coefficients of detennination (R ), root mean square errors (RMSE, g C m' d' ) and coefficients of variation (CV) 
for simulated GPP with the Vis using two options: original unsealed Vis versus scaled Vis, compared to tower daily GPP 
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